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I. INTRODUCTION 



The hyperfine structure of hydrogenic atoms is one of the well-understood problems both 
experimentally and theoretically. Especially, the muonium has played an important role in 
the precision test of QED because its radiative corrections have been calculated to high 
orders and its hyperfine splitting has been measured very precisely ||T|: 

Az/(exp) = 4 463 302.88 (16) kHz (0.036 ppm). (1) 

Furthermore, a new experiment is in progress to improve the precision of Az/(exp) to about 
0.007 ppm f2j. To match this experimental accuracy, it is necessary to improve the theory of 
the a 2 (Za) and a(Za) 2 non-recoil radiative corrections as well as the leading ln(Za) terms 
of order a 4 ~ n (Za) n , n = 1, 2, 3, and some relativistic corrections. This paper presents 
details of the calculation of the a 2 (Za) radiative correction. A preliminary report of this 
work has been published ||. 

As is well known, the bulk of the hyperfine splitting can be explained by the nonrela- 
tivistic quantum mechanics and is given by the Fermi formula [|J 
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where R^ is the Rydberg constant for infinite nuclear mass, and m and M are the electron 
and muon masses, respectively. 

Many correction terms have been calculated over several decades since the pioneering 
work of Fermi. Unfortunately, different terms were often evaluated by different methods 
making comparison of the results nontrivial in some cases. This causes a particularly difficult 
problem in identifying and evaluating higher order correction terms. Recently, however, 
Lepage and his collaborators have developed an approach, called NRQED, to deal with 
the nonrelativistic and weakly coupled bound systems consistently, starting from quantum 
electrodynamics (QED) 0-0. This provides a solid framework for evaluating higher order 
radiative corrections systematically and unambiguously. This series of papers deal with a 
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treatment of radiative corrections of the muonium hyperfine structure within the framework 
of NRQED. 

Before describing our calculation, let us summarize the previous results on the muonium 
hyperfine splitting Av. It is customary to classify the QED corrections to Av into three 
types: radiative non-recoil correction, pure recoil correction, and radiative-recoil correction. 
We use the convention such that electron charge is e and the charge of the positive muon 
is — Ze. Of course Z — 1 for the muon, but it is kept in the formula in order to identify 
the origin of corrections. Note that each radiative photon on the electron-line contributes 
a factor a, that on the muon line a factor Z 2 a, and one jumping from electron to muon 
a factor Za. This factor also arises from the effect of binding on the velocity distribution 
of atomic electrons. In addition, there are small corrections due to the hadronic vacuum 
polarization and weak interaction effects. Thus one may write 

Az/(theory) = Az/(rad) + Az/(recoil) + A z/(rad- recoil) 

+ Az/(hadron) + Az/(weak). (3) 

Purely radiative terms of orders a(Za) and a(Za) 2 have been known for some time ||: 

/ 3 5 
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Here a e and a M are the anomalous magnetic moments of the electron and muon, respectively. 
The appearance of the factor (1 + a M ) in (^) is in accord with our definition of Ep in 
Note that the number 14.88 in the a(Za) 2 correction is different from 15.39 reported in Ref. 
0. This is due to the recent discovery of two mistakes in the literature. The first error is 
in the calculation of the aiZa) 2 correction due to the vacuum polarization insertion in the 



transverse photon in Ref. || . Recently several people independently found ||iO|-4T2l that this 
contribution is E F a(Za) 2 4/5), not E F a(Za) 2 /n(— 2/3) given in Ref. 0. The second 
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error was caused by omission of a part of the contribution due to the vacuum polarization 
insertion in the Coulomb photon, Epa(Za) 2 /it(— 8/15) ln2, when it was combined with the 
radiative photon contribution EpaiZa) 1 /7r(15.10 ± 0.29) [ |TTf . 
The known recoil corrections add up to M 
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where 7 = Zam r , m r = mM/(m + M). The radiative-recoil contributions, which arise from 
both electron and muon lines and from vacuum polarizations, are given by [] 
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The a(Za)(m/M) and Z 2 a(Za)(m/M) terms are known exactly [§,|T! 
parts of the a 2 (Za) term were evaluated by Eides et al. flMf . 
The hadronic vacuum polarization contributes |l5j 



(6) 

The In 3 and In 2 



hadron = v - ; — =-(3.75 ± 0.24 )E F 
= 0.250 (16) kHz , 



(7) 



where is the charged pion mass. 

Finally there is a small contribution due to the Z° exchange. Our re-evaluation of the 
standard-model estimate |]TB|,|TT| gives Q 



x Eq. (6) of Ref. ||] is valid only for Z = 1 although it does not affect the muonium. We thank B. 
N. Taylor and P. Mohr for pointing out this oversight. 

2 This is in agreement with the corrected value given in Ref. [17| and has a sign opposite to that 
of Ref. [^]. The same result was also obtained by J. R. Sapirstein and by M. I. Eides. We thank 
B. N. Taylor and P. Mohr for calling a possible problem of sign to our attention. 



TABLE I. Contributions of various terms to the hyperfine splitting of the ground state muo- 
nium. (The new result of this paper is not included.) They are represented in units of kHz. The 
contribution from the muon anomalous magnetic moment is included in each non-recoil radiative 
correction term in the left column. 
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-0.065 



A v weak = —Gf tip 

~ -0.065 kHz. (8) 

Numerical values of terms given by Eqs. - (^[) are summarized in Table |I[ If one uses 
the value of a, Roo and M/m from Refs. |18| , |T9[] and |]]: 



a 



= 137.035 997 9 (32) (0.024 ppm) 



#00 = 10 973 731.568 30 (31) m" 1 , 

— = 206.768 259 (62), (9) 
m 

theoretical prediction for the hyperfine splitting of the ground state muonium, sum of the 
contributions listed in Table [J is given by 

Az/(old theory) = 4 463 302.27 (1.34) (0.21) (0.16) (1.00) (10) 

where the first and second errors reflect the uncertainties in the measurements of m M and 
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FIG. 1. Representative diagrams contributing to the a 2 (Za) radiative corrections to the 
muonium hyperfine structure in which two virtual photons are exchanged between e~ and The 
muon is represented by x . 



a^ 1 listed in (|9]). The third error is purely theoretical and dominated by the uncertainty in 
the last a(Za) 2 term of (^). The last one, about 1 kHz, is an estimated contribution from 
the order a 2 (Za) correction in A(rad). 

As is clear from ( |TUD one must know the a 2 (Za) pure radiative correction in order to 
improve the theoretical prediction further. Fig. [l] shows typical diagrams contributing to 
this order. Recently, terms represented by the diagrams (a) - (e) of Fig. [I] have been 
evaluated by Eides et al. ||20|| . Their results are as follows: 
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Az/(Fig.l(d)) = -0.310 742 •• • ^-^—^-Ep 



71 



-0.171 4 kHz, 



(14) 



where z = ln((l + \/5)/2). The results (|TTD, ([12]) and ([i~3| ) are analytic, while (|T4l ) was 
evaluated numerically after reducing the integral to one dimension. We confirmed these 
results by an independent numerical calculation. However, our purely numerical evaluation 
of Fig. 1(e): 



Ai/(Fig.l(e)) = -0.472 48 (9) 
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-0.260 6 kHz 



(15) 



disagreed with the semi-analytic result of Ref. [21]. With our help, Eides [22] found an error 



in the Table after Eq. (23) of Ref. 21]. Their corrected value is in good agreement with 



(0)- 

Fig. ^] shows the complete set of Feynman diagrams of the type (f) of Fig. [1|, which 
have not been evaluated before our work ||. The preliminary result of our calculation for 
all diagrams of Fig. El was 



Ai/(Fig.l(f)) = -0.63 (4) 



a 2 {Za) 



71 



-0.347 (22) kHz, 



(16) 



where the error is mainly due to the uncertainty in extrapolating the integral to zero infrared 
cutoff. The main purpose of this paper is to report a further improvement of this result: 
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FIG. 2. Two-photon exchange diagrams with fourth-order radiative corrections on the elec- 
tron-line. Diagrams which are related to these diagrams by time reversal are not shown explicitly. 
The muon is represented by x. 
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As a consequence of this result, the total contribution of the a 2 (Za) correction to the 
muonium hyperfine splitting becomes 

a 2 (Za) 

Aiz(Fig.l) = 0.767 9 (79) 1 Ep 

i\ 

= 0.423 5 (44) kHz. (18) 

This removes the dominant theoretical uncertainty in Ai/(theory). 

In Sec. m we outline the NRQED treatment of two-body bound system. It serves as 
the theoretical basis for the calculation of the a 2 (Za) correction as well as the calculation 
of the a(Za) 2 and higher order corrections discussed in the subsequent papers. In Sec. [I I 



we illustrate the general procedure of NRQED choosing the well-known aiZa) non- recoil 
radiative correction as an example. In Sec. [IV] we present our calculation of the a 2 (Za) 
purely radiative non-recoil correction to the muonium hyperfine structure. Some problems 
encountered in the numerical work are also discussed there. Sec. [V] is devoted to the 
discussion of our results. 



II. NRQED 
A. Why NRQED ? 

The Lorentz invariance has been one of the most important guiding principles for the 
development of quantum field theory. However, relativistic quantum field theory is often 
very cumbersome to apply to nonrelativistic bound systems. Such a calculation tends to 
be very complicated and requires an enormous effort, while the result reflects mostly the 
nonrelativistic feature of the system. For such a system an approach that incorporates most 
of the bound state effects from the beginning would minimize the amount of computation 
necessary to achieve the desired precision. For the case of electromagnetic interaction, this 
has been realized by a theory called nonrelativistic quantum electrodynamics, or NRQED. 
The NRQED enables us to avoid some, if not all, of the problems encountered in the usual 
treatment based on the Bethe-Salpeter equation. 
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The NRQED, formulated by Caswell and Lepage ||, is a rigorous adaptation of QED 
to bound systems. This theory enables us to take a consistent and systematic approach to 
loosely bound nonrelativistic systems. Compared with the conventional bound state theories, 
it allows easier power counting, more transparent cancellation of UV and IR divergences, 
and is manifestly gauge invariant. In spite of its superiority, however, the details of the 
theory has not yet been fully worked out. In this series of papers, we present an explicit 
construction of the NRQED Hamiltonian and develop a bound state perturbation theory 
based on it. 

As for the computation of the aiZa) and a 2 (Za) corrections, NRQED or any other 
relativistic bound state formalism gives the same simple recipe: calculate the forward scat- 
tering amplitude in QED and multiply it with |0(O)| 2 , where 0(0) is the nonrelativistic wave 
function at the origin. In NRQED this recipe can be directly justified by inspection of 
relevant diagrams and power counting. In other bound state formalisms, the corresponding 
procedure may be less straightforward. The latter approach becomes very complicated for 
higher-order corrections such as the a(Za) 2 and a(Za) 3 corrections. Difficulty in achieving 
high numerical precision by this method is one of the sources of theoretical uncertainty at 
present g|23|. 

The approach adopted by the NRQED, however, loses its effectiveness for the high Z 
system. In such a case it is desirable to avoid expanding in v ~ Za. Recently, an attempt 
has been made to calculate the order a term without expanding in Za p4}| . However, this 
approach may have difficulty in providing a good precision for Z = 1. This is primarily be- 
cause, for low Z systems, the bound electron is almost on-the-mass-shell and causes the near 
infrared divergence. As a consequence the convergence of numerical integration deteriorates 
as Z decreases. As is shown in the subsequent papers, the NRQED method enables us to 
deal with the near infrared divergence problem order by order in a systematic expansion in 
Za, and allows us to calculate the expansion coefficients with high precision. This is why 
the NRQED method is a powerful tool for low Z systems. 
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B. Outline of NRQED 



In the NRQED approach to the bound state problem, one first derives the NRQED 
Lagrangian from the QED, and then uses it to determine the correction to the energies and 
wave functions by a systematic application of the Rayleigh-Schrodinger perturbation theory. 

The NRQED Lagrangian consists of all possible local interactions satisfying the required 
symmetries, such as gauge invariance, parity invariance, time reversal, galileian invariance, 
hermiticity, and locality. We use the same photon Lagrangian (— 1/4) F^F^ as that of 
QED. In addition, new photon interaction terms are introduced to represent the insertion 
of the fermion loop, such as vacuum polarization and light-by-light scattering. 

In order to define the NRQED Lagrangian precisely, we must regularize the interaction 
terms of NRQED, e.g., by cutting off contributions of large momenta. Since this theory 
is meant to apply to nonrelativistic systems, the cut-off A may be chosen as the typical 
mass scale of the system, e.g., the rest mass of an electron. With the cut-off A thus fixed, 
the theory becomes well-defined, even though the interaction terms are strongly dependent 
on the cut-off parameter. In the following the cut-off is understood implicitly, and will be 
exhibited only when it is necessary. The choice of the momentum cut-off used for the NRQED 
scattering amplitudes is arbitrary but the physical quantity computed should be independent 
of any particular choice. In other words, the NRQED theory must have reparametrization 
invariance with respect to the choice of cut-off. This is analogous to the existence of the 
renormalization group in the renormalizable relativistic field theory. It is important to note 
that the NRQED is fully equivalent to the QED. The only difference is that it is better 
adapted to low energy bound systems. 

The NRQED rule for determining the operators which appear in its Lagrangian and their 
coefficients is simple and straightforward: Each term of the scattering amplitude calculated 
in the NRQED must coincide with the corresponding scattering amplitude of the original 
QED at some given momentum scale, e.g., at the threshold of the external on-shell particles. 
The center of mass frame is used for both bound state and scattering state calculations. Since 
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the same argument about reparametrization invariance holds for the momentum scale chosen 
for comparison of QED and NRQED scattering amplitudes, the at-threshold condition is just 
for convenience. However, the on-shell condition for the external fermion is more than a 
matter of convenience. In order to regulate the IR singularity it is convenient to introduce 
the photon mass A in the calculation of scattering amplitude of both QED and NRQED. This 
finite photon mass together with the on-shell condition ensures that the NRQED scattering 
theory has a pole in the region of the complex energy plane of the external fermion in which 
the scattering theory can be analytically continued to the off-shell bound state theory. 

We use the normalization ivu = 1 for the external 4-component spinors in the QED cal- 
culation instead of the conventional relativistic normalization uu = 1 so that both QED and 



NRQED S- matrix have the same normalization p5[ . This ensures that physical quantities, 
such as decay rate and cross section, calculated in both theories are the same. 

Note that the scattering amplitude of QED is fully renormalized, namely, it is finite and 
completely determined within QED. This enables us to fix the NRQED "renormalization" 
constants without ambiguity. This also means that the coupling constant a and fermion 
masses in the QED are the renormalized ones determined on-shell, and these a and fermion 
masses are used as the "bare" coupling constant and "bare" masses of NRQED. 

It is convenient to write the NRQED Lagrangian in two parts: L ma ; n and -^contact- The 
L mam part consists of the fermion bilinear operators. Fermions in NRQED are expressed by 
the Pauli two component spinor field ip(t,x) (instead of the Dirac spinor). If one takes into 
account the required symmetries of the theory, the main part of NRQED Lagrangian -L main 
must have the general form |].[(| 

u 3 2 3 4 

£main = V A + ~ ^ 



2m 8m 3 

ea-B e(D-E-E-D) 
+CF-: + c D - 



2m ' " 8m 2 
ieo ■ (3 x E - E x 3) e{3 2 ,a ■ B} 

+CS 5 h C W1 

8m z 8m A 
-eD l a-BD l e(a ■ DB ■ 3 + D ■ Bo ■ D) , , . . 

+C "-' 4m' + °>'' ^ + ' ' ^ ' (19) 
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where D t = d t + ieA° and D = d — ieA. (We put c = 1 and H — 1 henceforth.) The positron 
part can be written down in a similar way. The particle-antiparticle mixed interaction is 
not present in L ma ; n . The first three terms are related to the kinetic term of the QED 
Lagrangian. The second and third terms are derived from the expansion 

P p 



E = Jp + m 2 = m + ^-^—+ ... . (20) 

These three terms of (O) have coefficients unaffected by the radiative correction as a conse- 
quence of the renormalizability of QED, while the coefficients q of other terms are modified 
by the QED interaction and can be expressed as a power series in the coupling constants a 

Cj = c ; o) + c i {1) a + cf ) a 2 + ... . (21) 

Some of the operators in ((L9|) can be generated by the Foldy-Wouthuysen-Tani transfor- 
mation of the Dirac Lagrangian. These operators have the coefficient = 1 while other 
operators have cj ^ = 0. Note that c^'s do not have coefficients involving Za caused by the 
binding effect because they are determined solely by comparison of the NRQED and QED 
scattering amplitudes without referring to the bound states. 

Eq. (|TJJ) has an infinite number of terms. Not all of them, of course, are needed in 
a practical calculation. The operators necessary to carry out a particular calculation are 
determined by the power counting rule of NRQED for the bound state. We will show this 
process explicitly in the next subsections where the NRQED Hamiltonian is constructed. 

The L ma i n of NRQED alone is not sufficient to produce the same physical quantities as 
those from QED. To make NRQED equivalent to QED, we must add another term to the 
NRQED Lagrangian. It consists of terms of contact interaction type: 

^contact = dl ~M^^ ' (X^X) + ^(V^XX^) 

+ d 5 ^rf(^D 2 a^-(x^x)+ ••• , (22) 

where x represents a fermion field (of mass M) such as a muon or a positron. The third and 
fourth terms in fl22"|) are needed only when both electron and positron are present. This is 
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because, from the viewpoint of NRQED, the electron-positron annihilation is a high energy 
process and can only be represented as a contact interaction term. For the muonium, only 
the first and second terms are relevant. The fifth term is an example of contact terms 
including derivative interactions, which are of higher order in < p 2 /m 2 >~ (Za) 2 . 

The coefficients di are chosen such that these contact interactions make up the difference 
between the QED electron-muon scattering amplitude and the corresponding NRQED scat- 
tering amplitude derived from the Lagrangian L main . This procedure enables us to determine 
the coefficients di completely. 

As is clear from the above discussion, these NRQED "renormalization" constants q and 
di have the parameter dependence 



Of course, the experimentally observable result of calculation must be independent of the cut- 
off A, and gauge invariant. This is realized by a systematic application of the nonrelativistic 
Rayleigh-Schrodinger perturbation theory to the bound states. Note also that q and di are 
finite and well-defined in the infrared limit and hence require no infrared cut-off. 

Just as the actual execution of renormalization program of QED must rely on the covari- 
ant perturbation theory, a comprehensive formulation of NRQED can be realized explicitly 
only within the framework of the nonrelativistic Rayleigh-Schrodinger perturbation theory. 
This means that we have to choose an appropriate part of the Hamiltonian as the unper- 
turbed term and treat the rest as perturbation. 

To deal with the muonium we find it generally convenient to define the unperturbed 
system in terms of the ground state solution of the nonrelativistic Schrodinger equation: 



where m r is the reduced mass, E° = —^ 2 /{2m r ) is the ground state binding energy, 7 = 
[Za)m r being a typical momentum scale of the Coulomb bound state. The solution of this 
equation is 



q = Ci(a, A, to), 



di = di(a, Za, Z 2 a, A, to, M). 



(23) 




(24) 
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(j){p) 



(25) 



{p 2 + 7 2 ) 2 

The unperturbed electron field ip(p) is thus expressed by this wave function <p{p) times the 
Pauli spin factor. Using the remaining interaction terms in the NRQED Lagrangian together 
with the photon Lagrangian, we can construct the effective potentials. These potentials are 
to be treated as perturbation. 

The nonrelativistic Rayleigh-Schrodinger perturbation theory gives 



AE n = ifcV^l + H (Je V ) ^n]E= E o 
+ TPlV(G - 4^L)V*Pn \e=EO 



E ~ En' 



E-E2- 



E-E°- 



+ . . . . 



(26) 



tivistic 



The Green function Go{k, q : E®) appearing here is known in a closed form for the nonrela- 
tic Coulomb potential ||26|| . For the ground state n = 1, we find 

hm {Go-- r- Q —) = T -{27ry8 6 (k - q) 

—2m r —Ze 2 —2m r 



(£ 2 + 7 2 ) |£-g1 2 (<? 2 + 7 2 ) 



(27) 



where 



and 



R(k,q) = 



7 2 



2 fc 2 + 7 2 



-4- 



7' 



(fc 2 + 7 2 ) 2 (g 2 + 7 2 ) 2 
+ 5 l0g ^ + (4A-1)V 2 tan " 1(4A " 1)V2 



q 2 + 7 2 



A = 



(fc 2 + 7 2 )(g 2 + 7 2 ) 

A^ 2 \k — q\ 2 



(28) 



(29) 
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The first, second, and third terms of the expression (|27| ) can be understood as corresponding 
to zero, one, and two or more Coulomb-photon exchanges. 

In order to determine which terms of the Hamiltonian are needed to obtain the desired 
precision it is useful to know the expectation values of various operators with respect to 
appropriate wave functions PSI]. For a nonrelativistic Coulombic bound system, one finds 



< d >~ m(t>/c), < d t >~ m(v/c) 2 , < eA° >~ m(v /c) 2 , 

< eA >~ m(v/c) 3 , < eE >~ m 2 (v /c) 3 , < eB >~ m 2 (t>/c) 4 , (30) 

where m is the electron mass, v is the typical velocity of a bound electron, and c(= 1) is the 
velocity of light. Thus, in Eq. (|19D, the first two terms, next four terms, and the remain- 
ing terms correspond to the interactions which start at orders v 2 ,v A , and t> 6 , respectively. 
Radiative corrections, which alter the values of the coefficients q's and d^s, will keep the 
estimate (p0|) intact. The information (|30|) can be used to terminate the series of interaction 
terms at the desired precision. In this sense, the NRQED Lagrangian is an expansion in 
both the coupling constant a and the velocity v. 

The NRQED "renormalization" coefficients play important roles in restoring gauge in- 
variance which might have been broken by regularization. The explicit form of these coeffi- 
cients depends on the regularization method. Gauge invariant regularization is desirable but 
not necessary. If one proceeds carefully, even a simple momentum cut-off method may be 
used [p8| . (This is only true for an Abelian gauge theory such as NRQED.) In a calculation 
of the would-be divergent quantity in NRQED, we put the UV cut-off A not in the fermion 
momentum but in the photon momentum |29| . 



Because of the way NRQED is constructed, gauge invariance of the NRQED amplitude 
with its complete set of "renormalization" constants is automatically guaranteed by the 
gauge invariance of the corresponding QED amplitude. Since QED and NRQED are sepa- 
rately gauge invariant, we may choose different gauges for QED and NRQED. We will use 
the Feynman gauge for QED calculation, and we use the Coulomb gauge for NRQED. The 
Feynman gauge minimizes the amount of work for numerical computation, and the Coulomb 
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gauge is more suitable for describing the nonrelativistic behavior of the electron. 

The dominant contribution to the hyperfine splitting between the spin J=l and J=0 
states originates from the interaction between the electron spin and the muon spin mediated 
by a transverse photon of momentum k. This leads to a potential of the form 

Vf = 2^ * x m ' w (xt(_fc) x (31) 

in the momentum space representation, where M is the muon mass. Using this Fermi 
potential Vf in the first order perturbation theory and taking the difference between J=l 
and J=0 states, we obtain the hyperfine Fermi splitting 

E F = (n = l\V F \n = l)\j=l 

d 3 p f d 3 k {8^¥f) 2 



(27r) 3 J (27r) 3 (p2 +7 2)2(|^ + ^|2 +7 2 ; 

-^—({k x <j e ) ■ ^ 6 ((— k) x a„))^— 
2m u ; 2M u ; m;/ k 2 

8(Za)7 



2 



3 



3 mM 

Needless to say, the Fermi potential is of order v A (m/M)m ~ (Za) 4 (m/M)m. 



(32) 



Various interaction terms and propagators are represented by the NRQED "Feynman" 
diagrams shown in Fig. [|. It is convenient and useful to express higher-order amplitudes 
representing scattering states or bound states by corresponding diagrams. 

We classify the diagrams according to the number of external photons in the QED Feyn- 
man diagrams. In the following subsections we shall show step by step how the corresponding 
NRQED Lagrangian L main (or Hamiltonian H ma ^ n ) is determined. 



C. Scattering by a static external potential 

We have already found the general form of the main part of the NRQED Lagrangian 
given by Eq. ( |I9D using the required symmetries for the theory and the power counting rules. 
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Propagators 

Coulomb photon propagator 



q 2 + X 2 



Transverse photon propagator 



q 2 +\ 2 



(q ) 2 - q 2 - X 2 + it 



Fermion propagator 



1 



Vertices 



2m 



Coulomb vertex 



Dipole vertex 

p +p 

—e 

2m 

A • A vertex 

e 2 S ij 



2m 



FIG. 3. NRQED "Feynman" rules for vertices and propagators. They can be used for both 
scattering and bound state calculations. E in the fermion propagator represents the fermion's 
bound state energy: E = for scattering and E = — j 2 /(2m r ) for the ground state muonium. The 
photon mass A is set to zero in the bound state calculation. 
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Figure 3 (continued) 



Darwin vertex 



8m 2 
Fermi vertex 



Relativistic kinetic vertex 



Seagull vertex 
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Spin-orbit vertex 
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Time derivative vertex 
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Figure 3 (continued) 

Derivative Fermi vertices 



p'p vertex 



cvp 



m? 



tvp 



Therefore, the remaining task for construction of the NRQED Hamiltonian is determination 
of the coefficients of these operators appearing in Eq. (p!9|). 

Let us first consider the QED diagram in which one photon is exchanged between the 
electron and the muon. The first step to obtain the "renormalization" coefficients of the 
operators in the NRQED Hamiltonian is to carry out nonrelativistic reduction of the QED 
scattering amplitude exchanging one photon between the electron and the muon. Compar- 
ing this QED scattering amplitude with the scattering amplitude derived from the general 
form of NRQED Lagrangian given in Eq. (0), we are able to fix the "renormalization" 
coefficient Cj's. We want to chose the simplest process to find them. It turns out that all 
"renormalization" coefficients in can be obtained by using the external static potential. 
The comparison between the corresponding QED and NRQED amplitudes is shown in Fig. 

20 

1 



|j. We will work out nonrelativistic reduction of operators of order up to av 6 and av 4 for 
spin-flip and spin-non-flip ones, respectively. Since the spin-non-flip operators contribute 
to the hyperfine splitting only through the higher order bound state perturbation, we need 
only operators of order lower than the spin-flip ones. The QED scattering amplitude to be 
studied here consists of a tree vertex and one dressed by a radiative photon. Because we are 
dealing with the scattering amplitude, we also have a diagram with self-energy insertion on 
the external fermion lines. However, these diagrams can be dropped after the mass renor- 
malization and wave function renormalization are carried out, if one chooses the on-shell 
renormalization scheme. Then the QED scattering amplitude is expressed by the usual form 
factors Fx and F 2 . 

For the external static vector potential A(q), we easily find the QED scattering amplitude 



eu(p') 

= fi(?V(p' 



-7 ■ A(m(<l 2 ) + A\q)q^ F 2 {q 2 ) 

2m 



u(p) 



C -* -*■ 



{V + p 2 )a- (qx A) + 



8m 3 
+F 2 (q 2 )^(p' 

ie 



ip(p) 



ie 



ie 



(p +p )a-(qx A) 



8m 3 



a ■ p'a ■ (q x A)a ■ p + 



(33) 



where u and ip are Dirac and Pauli spinors, respectively. Similarly, for the external static 
Coulomb field A°(q), we have 



eu(p') 
= F 1 (q 2 )^(p' 
+F 2 {q 2 )^{p>) 



^A^qlF^q 2 ) - ^-a°JA°(q)q>F 2 (q 2 ) 
2m 



u(p) 



eA° 



-A 



2 A° + -& 2 t-(p'xp)A° + ... 



ie 



Am 2 



ip(p) 



Am z 2m z 



ip(p) . 



(34) 



Taking account of the fact that q° is of order v 2 and |g| is of order v, the nonrelativistic 
expansion of the form factors can be written as [B7[ 



m 2 ) 

F2{q 2 ) 



a 
3^ 



m c 
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+ 0(av\a 2 v 2 " 



a q- 



+ 0{av\a 2 v 2 ) 



tc 12m 2 



(35) 
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FIG. 4. QED and NRQED scattering diagram comparison. The diagrams on the left and 
right of the = sign represents QED and NRQED diagrams, respectively. The external fermions are 
on-the-mass-shell and at threshold. Self-energy diagrams coming from the A ■ Atp^tp vertex as well 
as self-mass counterterms are not shown explicitly. 
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where a e = ^(0) is the anomalous magnetic moment of the electron. 
Combining (|33|) 



and fl35| ) together, and comparing with the scattering amplitude 
derived from the NRQED Lagrangian given by (|19D, we find that the "renormalization" 
coefficients must be chosen as 
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(36) 



This procedure enables us to construct the NRQED Hamiltonian. However, it does not 
provide the complete NRQED Hamiltonian. One must also include terms which are the 
NRQED analogues of QED counterterms such as —Smuu. To see this let us calculate the 
NRQED scattering amplitude which arises from the Coulomb term ifj^eA ^ modified by a 
1-loop NRQED radiative correction by a perturbative treatment of H ma i n . 

The perturbation here means that the zeroth-order of NRQED Hamiltonian contains only 
the free part of the electron, and thus the Coulomb interaction is treated as perturbation. 
Some of these scattering amplitudes involving radiative corrections require new forms of 
the NRQED operators while others may be represented by additional "renormalization" 
constants of the already existing operators in -ff ma i n - 

The fermion kinetic energy term in (|19D gives the interaction term —^e(p'+p)-A/(2m)ip 
in the NRQED Hamiltonian. Note that, although the NRQED Hamiltonian is not an expan- 
sion into multipoles, we call this term the dipole interaction in the following for convenience's 
sake. Thus we consider the Coulomb term dressed by the transverse photon with the dipole 
couplings. The NRQED Feynman rule applied to this diagram gives 



t , e_y. f A d*k 1 

\m) l J (2vr) 4 (A;0)2_£2_ A 2 + , e 



p ■ p 



p-kp'-k 

k 2 + X 2 
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-eA°- 



1 



(37) 



E + k° - (p> + kf/(2m) + ie E + k° — (p + k) 2 /(2m) + ie 
We chose the contour in the upper half k° plane to pick up only the negative energy photon 
pole. Then we neglect k in the kinetic energy term (p+k) 2 / (2m) in the electron propagators. 
This is justified because the energy transfer between electrons is of order v 2 when \p\ is of 
order v, while the space component of the photon momentum \k\ is of order v 2 . After this 
approximation, angular integration over the photon momentum k becomes trivial, leaving 
only the \k\ integration: 
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(38) 



In the last step, we used the on-shell, at-threshold condition, E = —p 2 /(2m) + 0(v 4 
The diagram with a self-energy on the external electron line gives 
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(39) 



Note that the term — \j\/k 2 + A 2 is the "mass renormalization term" of NRQED. f] 

In order to maintain the equivalence of QED and NRQED we must include the neg- 
ative of these contributions in H ma[n . (See Fig. U(d).) ^From ([38l) and (|39|) we see 
that this is achieved by adding the new "renormalization" coefficients to the Darwin term 
-^eq 2 A°/(8m 2 )ip: 



5BI' 



(40) 



3 "Tadpole" diagram due to the < A ■ Aip^ip > completely vanishes after mass renormalization 
because this diagram does not depend on the external fermion momentum. 
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where 



NRQED 
Z D 



a 8 

7T3 



I 2A / 6 



(41) 



The entire coefficient of the Darwin term is the sum of QED and NRQED contributions: 
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(42) 



Actually, this additional contribution from NRQED serves to eliminate the contribution of 
the longitudinal polarization associated with the finite photon mass ||30|| . In other words, 
the In A term in the "renormalization" coefficients due to QED is effectively replaced in the 
NRQED "renormalization" constant by 

5 



In A ->ln(2A) - 



6 



(43) 



Similarly the NRQED radiative correction to the Fermi term, —ipHea ■ (q x A)/ (2m)ip, 
yields the correct "renormalization" coefficient of the W\ and W2 derivative Fermi terms, 
ijjHe(p' 2 + p 2 )^ ■ (q x A)/(8m 3 )if} and ijjUe(—2p' ■ p)a ■ (q x A)/(8m 3 )ijj, respectively, which 
are given by 
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(44) 



The radiative correction comes also from vacuum polarization. Since vacuum polarization 
is a highly virtual process within the framework of NRQED, no vacuum polarization term 
exists in H main . Instead, its contribution is represented by the new photon interaction terms 
in NRQED. Again we begin with the nonrelativistic reduction of QED amplitude with one 
vacuum polarization insertion. QED gives the renormalized vacuum polarization tensor 



n'-'(g) = (g*V - sTfmq 



(45) 



with 
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u ^ = -&fo \-JX-e)-> - (46) 

For the second order, the photon spectral function p 2 (t) is known to be 

/N at 2 (l-\t 2 ) 

Expanding IT(g 2 ) around q 2 = 0, we obtain 

-»2 

n 2 (g 2 ) = c vp ^- + 0(av 4 , a 2 v 2 ) , (48) 

where 

*- = m ■ m 

Thus, in the Coulomb gauge, two new photon interaction terms are added to the photon 
Hamiltonian 

Cvp A\q)^(q)(^-^), (50) 
m A q A 

and 

c vp A°(q)^A°(q). (51) 



D. Photon- Fermion Scattering Amplitude 

Let us now turn to the processes which contain two fermion operators and two external 
photons. To determine the " renormalization" coefficients we must carry out the nonrela- 
tivistic reduction of these QED scattering amplitudes. In practice, however, we don't have 
to do it at all because the "renormalization" coefficients of these operators up to v 6 for spin- 
flip ones and v 4 for spin-non-flip ones are identical with those determined by the scattering 
amplitude due to a static external potential because of gauge invariance. For instance, the 
same "renormalization" coefficient Cs for the spin-orbit interaction term 

^^a.(p'xp)A°i; (52) 
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must be used for both the seagull term 



— i —* 

^t— ?.( g lxi( gi ))i°(^ (53) 
and the time derivative term 

(54) 

The seagull term is the only operator involving two-photons relevant to our immediate 
interest. This contributes to the {Za) 2 and a(Za) 2 corrections. 

In general an explicit nonrelativistic reduction of the photon-fermion scattering ampli- 
tude is necessary only if one wants to find the "renormalization" coefficients of operators of 
higher order in v, such as ip^E ■ Eip/m 3 ~ v 6 . 

We note that i? ma in is not an unique expression. Using the equation of motion for the 
fermion field, we can obtain another form of Hamiltonian. When -ff ma m is quantized, we 
should be more careful. Use of the equation of motion is equivalent to the transformation 
of the electron field. We have to take into account the Jacobian of this change of variables. 
Once the Jacobian is taken into account, two Hamiltonians become completely identical and 



produce the same results even for the bound state calculation |2"5(] . This is why we excluded 
the operators having time derivatives, such as ip(iD t ) 2 ip/m, from our consideration, since 
the equation of motion renders iD t ip to (p 2 /(2m) + (9(t> 4 ))■?/;. 

In this manner we have obtained all operators in the main part of the NRQED Hamil- 
tonian f/main necessary for our calculation to the desired order. 

E. The NRQED Hamiltonian F main 

For later reference let us write down the part of the NRQED Hamiltonian if ma m valid to 
order a by putting together the results of Sec. [11 C and [11 D| . To do this, we introduce the 



q 2 derivative Fermi term by combining the W\ and W2 derivative Fermi term at the order 
a. It is of the form: 
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(55) 



where p' and p are the outgoing and incoming fermion momenta, respectively, and q = (q°,q) 
is the incoming photon momentum. In the seagull vertex, q{ is the incoming momentum 
of the vector potential A. The superscript A indicates that the Hamiltonian is regularized 
with the UV cut-off A. The "renormalization" coefficients are 
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(56) 

Thus far we have not shown explicitly the contact term H contSuCt of the NRQED Hamil- 
tonian, which is also obtained by comparison of the electron-muon scattering amplitudes in 
QED and NRQED. The explicit form of the contact term will be given in Sec. Ill and IV 
as we calculate a(Za) and a 2 (Za) corrections, respectively, to the hyperfine splitting. 
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F. Application of H ma i n to Bound States 

Let us now turn our attention to the bound state calculation using H mSLin . The main part 
of the NRQED Hamiltonian for the muon field is obtained by replacing the charge e by — Ze 
in the H ma \ n for the electron field. For the nonrecoil hyperfine correction, only the Fermi and 
Coulomb terms are necessary in the muon Hamiltonian. Together with the photon Hamil- 
tonian, we can construct various perturbative potentials appearing in the nonrelativistic 
Rayleigh-Schrodinger perturbation theory (|2BJ). The lowest order contribution Ep to hyper- 
fine splitting comes from the Fermi potential (|HT|). A survey of Eqs. Q55D and ( |5ED shows 
that the only order a correction is a e Ep which exhibits the effect of the "renormalization" : 
Cf — 1 = a e - Other possible contributions to the hyperfine splitting coming from i? ma i n are 
those of the first order perturbation of the derivative Fermi term and the seagull term, and 
the second order perturbation which involves the Fermi term and the p A relativistic kinetic 
term or the Darwin term. The (p' 2 + p 2 ) derivative Fermi term leads to the potential of 
order (Z a) 6 (m / M)m: 

v - - --i^f 1 ^ ■ <^ x wwhv ■ (57) 

The Darwin term generates the potential of order (Za) A m: 

Vd = ^ Wx) ^_. (58) 

Expectation values of these potentials with respect to the bound state wave function 
diverge due to integration over q. This is why we need the help of the contact term H contact 
for their cancellation. When the effect of H contact is included, these four potentials together 
give the {Za) 2 Breit correction. A detailed discussion about the treatment of these UV 
divergent operators is found in [10 where the derivation of the Breit {Za) 2 correction from 
NRQED is described. 

Similarly an a(Za) 2 correction is obtained when the contribution of the 
"renormalization" coefficients is included in each potential. Third order perturbation theory 
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in -f/main with an intermediate radiative photon and dipole couplings also gives the a(Za) 2 
correction. This is because these diagrams have the structure similar to the derivative 
Fermi term or the Darwin term as is shown in the determination of the "renormalization" 
coefficients in NRQED (See Eqs. fl38|) and (0)). 

The additional photon interaction terms ( |5D| ) and fl5I|) due to vacuum polarization pro- 
duce the effective potentials 

Vtvp = ^M~^h^^ ^ ' iX ^ X ^ X \q2 q +\2)2 ( 59 ) 

and 

^v P = "^^Wx) • (60) 

We note that the first spin-flip potential V tvp has exactly the same structure as the q 2 
derivative Fermi potential. Thus it contributes not to the order a(Za) but to the order 
a(Za) 2 . The spin-non-flip potential V cvp behaves as a S function potential in the coordinate 
space just like the Darwin potential. Thus it also contributes to the a(Za) 2 term through 
the second order perturbation theory. 

To summarize, no order a correction exists besides a e Ep, where Ep is given by fl3"2|). In 
the NRQED formulation, it is transparent why only the anomalous magnetic moment of a 
free electron contributes to the order a correction to Ep. 

III. THE a(Za) CORRECTION 

In this section we show how the non-recoil radiative correction of order a(Za), calculated 
long ago by Kroll and Pollock and by Karplus, Klein and Schwinger | 3lf , can be obtained 



within the framework of NRQED. For brevity, let us refer to this as the K-P term. The 
procedures developed here are readily applicable to the a 2 (Za) term calculation in NRQED. 
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(a) (b) (c) 

FIG. 5. QED diagrams contributing to the a(Za) radiative correction to the muonium hy- 
perfine splitting. 



A. Diagram Selection 

The QED diagrams involved in this calculation are shown in Fig. [5[ In the original 



and subsequent works [^T|,^,^,|T^] , the K-P a{Za)Ep pure radiative correction was evalu- 
ated from the QED diagrams with the external fermions put on the mass-shell and at the 
threshold, and multiplied by the square of the nonrelativistic Coulomb wave function at the 
origin. This recipe was justified after complicated and rigorous consideration of the rela- 
tivistic bound state theory. We shall show that NRQED provides an alternative justification 
of this procedure in the sense that no other correction term is needed in this order. 

The correction terms whose coefficients are odd powers of Za may arise only from very 
limited sources in the NRQED bound state theory. The NRQED Lagrangian L ma m consists 
only of terms with even parity. This implies that the expectation values of these terms with 
respect to the Coulomb wave function are even in Za, the typical electron momentum of 
the Coulomb bound state being \p\ ~ [Za)m. 

The odd power of Za in the K-P term ~ a(Za) 5 m 2 /M therefore implies that there 
is no contribution to it from the L ma m part of the NRQED Lagrangian. [Note that the 
"renormalization" constant Cj in L main does not depend on Za. (See Eq. ( |T9| ) and Eq. 

(ED)-] 

This means that the correction we are looking for must come entirely from the NRQED 
contact terms of (^2|) . To determine the contact term, we compare the scattering amplitudes 
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FIG. 6. QED and NRQED two-photon exchange scattering diagram comparison in the pres- 
ence of the radiative correction. The shaded circle represents the contact term introduced in 
this comparison. The NRQED diagrams in the bottom lines actually contribute to the a(Za) 2 
correction. 



evaluated in NRQED and QED in the same power of explicit a and Za. This comparison 
is shown in Fig. |^. For a given power of the coupling constant a, the number of QED 
diagrams is finite while the number of NRQED diagrams is infinite. We terminate the series 
of NRQED scattering diagrams using the power counting rule for their contribution to the 
bound state. 

We chose the electron mass m as the momentum scale of comparison, and evaluate the 
scattering amplitudes of both QED and NRQED on the mass-shell and at the threshold. In 
general, this procedure must be carried out for both spin-flipping and non-flipping ampli- 
tudes. However, for the K-P term, only the spin-flipping one is needed. The spin-non-flipping 
type produces a Lamb-shift type contact term, which contributes to the hyperfine splitting 
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only in the order a(Za) 7 m 2 /M and above. 

As we have discussed in Sec. [TI], the comparison of QED and NRQED scattering ampli- 
tudes gives rise to a contact term to the NRQED Hamiltonian. We restrict ourselves to the 
consideration of the contact term relevant to the hyperfme splitting, i.e., 

5H = -^-Lty+o^) . (x V- x); (61) 

because this is the only source of the K-P term as was discussed above. 

Let us first focus on the contribution from the vacuum polarization insertion. The two- 



photon exchange scattering amplitudes containing the vacuum polarization potential of fl55|) 
contributes to the order a(Za) 2 , not to a(Za). Thus the only contribution from the vacuum 
polarization is obtained from the contact term which is determined by calculating the QED 
two-photon exchange amplitude with one vacuum polarization insertion in the photon line 
with the on-shell at-threshold external fermions times the square of the Coulomb wave 
function at the origin. 

Let us turn next to the contribution from the radiative photon. The QED diagrams 
related to this correction are shown in Fig. [5[ All three QED scattering amplitudes have 
the same form: 

■rrQED _ 2(7 2\2 f d 4 q U e £^ v U e U m M.^ v U m 

11 ~ e [Z€ > J (2vr) 4 {q 2 + ie) 2 • [b ' 

Here the electron factor £^ v is different for each diagram but the muon factor M.^ v is common 
to all these diagrams and represents the sum of the ladder and crossed-ladder diagrams: 

jjJJ- A + M) lv 7 „(/+ ^ + M) 7m 
u (I - qf -M 2 + ie (1 + q) 2 — M 2 + ie ' 1 ' 

where I = (M, 0) is the external muon momentum and q is the four momentum flowing in 
the loop between the electron and the muon. As is well known ||, in the limit of infinite 
muon mass, the muon factor reduces to 

AV = 7 , Al» 2M W 1 ■ (64) 
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FIG. 7. Ladder and crossed-ladder diagrams. 



This is represented by the symbol x in Figs. [I], ^ and |. 

The hyperfine splitting projection operator for strings of 7 matrices is obtained by taking 
the difference between the J = 1, J z = state and J = state and using the spherical 
symmetry of the system fl3"3"[ : 

1 ^Tr^W + IF^AVtW + 1)] , (65) 
iZ i=i 

where Roman letters run from one to three while Greek letters run from zero to three. This 
projection is true only for external fermions on-the-mass-shell and at-threshold. The trace 
of the muon factor is easily taken, yielding 

, -2ni5(q ) 

e ^ iq 2M ■ (66) 

We take the e^ji g- 7 part together with the electron projection operator as the hfs projection 
operator, and the other muon factor will be included as a numerical factor. 

In order that these diagrams contribute to the hyperfine splitting, one of the exchanged 
photons must be transverse (attached to a vertex 7*) while the other is Coulombic (attached 
to a vertex 7 ). Our projection operator of hyperfine splitting picks up automatically this 
structure from the electron-line. 

The corresponding NRQED scattering amplitude consists of many diagrams, but most of 
them actually contribute to the order higher than the K-P term. The only diagram necessary 
is a combination of the Fermi potential multiplied by the NRQED renormalization constant, 
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namely the second order anomalous magnetic moment, and the Coulomb potential. This 
scattering amplitude is named iT NR Q ED . 

Other diagrams, such as the combination of the Darwin potential Vd including the 
"renormalization" constant and the Coulomb potential, have the same power of explicit 
a and Za as the Fermi one, but diverges linearly in both UV and IR region. These diver- 
gences cancel out in the bound state calculation. The detail is similar to the discussion on 



the Breit term calculation given in [K| . Eventually they contribute to the terms of order 
a(Za) 2 . This argument holds also for the potentials Vtvp and V cvp representing the vacuum 
polarization effect. 

The QED processes with three or more photon exchange contribute to obviously higher 
order terms due to the explicit extra power of the coupling constant Za. 

The contact term can be defined as the QED amplitude minus the NRQED amplitude 
for the two photon exchange process: 

- d^^a^) ■ ( X %x) = iT® ED - iT NR ® ED . (67) 

Actually both QED and NRQED amplitudes are IR divergent in the limit of the vanish- 
ing external photon momentum q. These threshold singularities cancel each other in the 
difference 

This contact term is to be put into the first order perturbation theory. Then the wave 
function integral is trivially done, resulting in the square of the Coulomb wave function at 
the origin. Thus the K-P term is given by 

Au(KP) = \m\ 2 ^(°e ■ aX=o , (68) 

where |0(O)| 2 = 7 3 /7r for the ground state. In the actual calculation, we take the difference 
between spin J=l and J=0 for the scattering amplitude first using the projection operator. 
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B. Calculation of the QED Amplitude 



We have shown that the a(Za) non-recoil radiative correction comes entirely from the 
NRQED contact term evaluated at the origin of the wave function. On the other hand, 
evaluation of the NRQED contact term is equivalent to that of the on-shell at-threshold 
QED scattering amplitude. This is why the calculation of the a(Za) term is much simpler 
than other terms such as the a(Za) 2 term. 

Our approach to carry out the computation of the QED scattering amplitudes is by nu- 
merical integration. Let us explain the outline of our procedure. The detail of the calculation 
is given in Appendix ||. The electron-line structure of each diagram is directly written down 
using the parametric Feynman-Dyson rules for QED [04 , 3q1 . Feynman parameters assigned 
to the electron-line are zi, Z2, and Z3, while one assigned to the radiative photon line is 
Z4. The momenta flowing in the fermion lines after the radiative photon loop momentum 
is integrated out are expressed in terms of correlation functions Bij, which are functions of 
Feynman parameters and determined by the topology of the loop structure of the diagram 
alone. Then our integrals are expressed as two or three dimensional Feynman-parametric in- 
tegrals with an additional one dimension corresponding to the magnitude of the momentum 
q of the external potential. 

Two of the QED diagrams have UV divergences and must be renormalized. The renor- 
malization terms are generated using the projection operators in the algebraic program 
FORM [p6| . Our projection operators for QED renormalization constants are quite general 
and applicable to any order. They are presented in Appendix [A|. All of the renormalization 
constants are determined in the on-shell scheme. These renormalization terms should be 
expressed by the same Feynman parameters as those assigned to the original diagrams in 
order to realize point-by-point subtraction in numerical integration by means of the adaptive 



iterative Monte-Carlo integration routine VEGAS [|37|j . 

The hfs contribution due to the second order anomalous magnetic moment should be 
subtracted from the diagrams involving the second order vertex correction. Actually it is 
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very easily done along with the charge renormalization: Let the external photon momentum 
q tend to zero in the original diagram expression of the electron factor. Then subtract this 
IR limit from the original diagram. We can easily prove that this IR limit of the diagram 
is nothing but the sum of the charge renormalization constant and the anomalous magnetic 
moment of the second order. (See [ID] for details.) 



Even though all diagrams are free from UV divergences after the renormalization is 
completed, they still suffer from IR divergence. In general, the Coulomb bound state has 
two kinds of IR divergence: one is due to the threshold singularity, and the other is due to 
the radiative photon. 

The mechanism of threshold singularity is the following. In order to contribute to the 
hyperfine splitting, one of the two exchanged photons must be Coulomb-like while the other 
is transverse. This Coulomb photon may be absorbed in the wave function. As a result, 
the diagram is reduced to one of lower order in Za, or multiplied by l/(Za). This is 
the physical origin of this type of IR divergence, which is ubiquitous in the relativistic 
treatment of bound state problem. In the calculation of the a(Za) correction, however, 
such a "divergence" can be avoided completely by subtracting the contribution of the free 
anomalous magnetic moment. This is because other threshold singularities are absent due 
to the on-shell renormalization. 

The remaining IR singularity is caused by radiative photons. Our choice to deal with this 
singularity is to put a small photon mass A in the radiative photon. For the K-P term, this 
singularity must cancel out when all QED diagrams of the gauge invariant set are included. 

C. Summary of the a{Za) Correction 

We have shown that non-recoil radiative corrections to the muonium hyperfine splitting 
having the odd power of Za comes only from the contact term of NRQED, and that this 
contact term is determined as the difference between the QED and NRQED scattering 
amplitudes. 
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The resulting expression for the K-P radiative correction can be evaluated either ana- 
lytically or numerically. We have chosen the later approach. The three dimensional inte- 
gration has been carried out using the adaptive iterative Monte-Carlo integration routine 



VEGAS [0. Each diagram has the IR divergence of the form proportional to ym/X, but 
their sum is finite. Our numerical evaluation shows that the contribution due to the radiative 
photon is given by 

Az/(KP) ph = -2.556 80(6)a(Za)E F . (69) 

We have also evaluated this integral analytically and obtained the same result as that of 
Kroll and Pollock, and Karplus, Klein, and Schwinger [3j| : 



Az/(KP) ph = (In 2- — )a{Za)E F = (-2.556 852 • • -)a{Za)E F . (70) 



4 J 

An easy analytic calculation of the vacuum-polarization contribution gives 

Az/(KP) V p = -a(Za)E F . (71) 

Putting these results together we obtain the well known aiZa) correction in the frame work 
of NRQED: 

Az/(KP) = fin 2 - — + - J a{Za)E F . (72) 



This justifies the procedure adopted in Ref. [ 3 1 . 



IV. THE a 2 {Za) CORRECTION 

A. Diagram Selection 

In this section, we give an outline of the evaluation of the a 2 (Za) correction to the Fermi 
frequency E F which comes from the six gauge invariant sets of QED Feynman diagrams 
represented by Fig. [TJ. Our treatment of the bound state to find the contribution to hyperfine 
splitting coming from these diagrams is completely identical with that of the a(Za) K-P 

38 



correction. A new diagram appearing in this order is the light-by-light scattering insertion. 
The light-by-light scattering is a high energy process in NRQED. Thus it is represented only 
by a contact term in NRQED. As a result we have to include the four-photon interaction 
in the NRQED Hamiltonian. But as an operator it contributes to orders higher than our 
interest here. Therefore, what to do is again to calculate the contact term starting from the 
scattering amplitudes of these diagrams with the on-shell at-threshold particles and then 
subtract the contribution of the fourth order anomalous magnetic moment from Figs. |l](d) 
and |(f). 

The numerical evaluation of Fig. [I] (a) - (e) can be carried out easily and our results 
are consistent with those previously obtained by Eides and his collaborators p0|-|22|j. In 
contrast, the diagrams of Fig. [I] (f) require a substantial effort to compute. A complete 
evaluation of this contribution is the main result of this paper. 

B. Calculation of the QED Amplitude 

Let us now discuss some technical details of calculation of (|T7| ) represented by the nine- 
teen diagrams of Fig. |2|. Since the bound state structure of these diagrams is identical with 
that of the a(Za) correction, the procedure of numerical evaluation of the a(Za) correction 
given in Appendix |B| can be applied readily to these diagrams. We applied numerous tech- 
niques developed for the numerical calculation of the anomalous magnetic moment g — 2 of 
the electron ||35|| , except that we avoided the use of "intermediate" renormalization which 
was introduced in the g — 2 calculation to avoid the IR singularity of each diagram. Instead 
we use the conventional renormalization procedure which is IR singular in the radiative 
photon mass A. This is because these IR divergent terms are needed to cancel out the other 
IR singularity, the threshold singularity in the vanishing external photon momentum q = 0, 
in the proper diagram. The detail of this mechanism is described in the calculation of the 
a(Za) K-P correction in Appendix [B|. 

It is convenient to divide the nineteen diagrams into four groups: 
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Group 1: Diagrams containing fourth-order vertex corrections. They are represented by the 
diagrams H 01 , H 02 , H 03 , H 09 , H w and H u of Fig. |. 

Group 2: Diagrams containing fourth-order self-energy insertions. They are represented by 
the diagrams if 04, H\ 2 of Fig. |2|. 

Group 3: Diagrams in which radiative photons span over two external photons. They are 
represented by the diagrams H 05 , H oe , H 07 , H 08 , Hi 3 , H u , Hi 5 , and H w of Fig. |[. 
Group 4: Diagrams containing two non-overlapping second-order radiative corrections. They 
are represented by the diagrams if 17, His, and if 19 of Fig. |[ 

The integrands corresponding to the individual diagrams of Fig. ||] were initially gen- 
erated using the algebraic program SCHOONSCHIP f38|. Later we generated the same 



integrands by FORM [36 check. 

The parametric representation of Group 1 diagrams is of the form 



3 a 2 (Za) m y°° dq 

32 7T 7T JO <J 



{dz) 



1-4 



1+ jq + l 



r 



u 2 V 2 

where the diagram Hqi, for example, has the electron-line operator 



(73) 



(74) 



Other diagrams of this group are obtained by permutation of 7 matrices. (See Ref. ]3S] for 
the definition of U, V, Di, etc.) Using the hyperfine splitting projection operator, one finds 
that the terms contributing to the hyperfine splitting are proportional to at least q 2 , and 
kills one of the <f 2 's in the denominator in Eq. ([73]). Thus Eq. ([73]) leads to the energy shift 
of the form 

dq f (c?2;) 1 „4 



Az/ G1 = — ^ ! -E F — 



V 2 uv U 2 V° 



(75) 



7T TT" JO {-q")J U 2 

where 1/V° is a symbolical representation of — In V in which the UV divergence is regularized 
and subtracted by the corresponding counterterm. (See Appendix [FJfor a precise definition.) 
The parametric representation of Group 2 diagrams is of the form 



3 a 2 (Za) 

32 7T 



m 



dq 



E F - / ^ 

n z jo q 2 



H- a + m 



{dz) 



2-4 



U 2 V 



7+ jq + m 



-q 



7 



(76) 
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where, for example, the diagram i?04 has the electron-line operator 

F 2 = 7 "( A + m) 7 /3 ( A + m)j ( A + m) la • (77) 
Its contribution to the hyperfme splitting has the form 

, ,rF„ W, i 

(78) 



a 2 (Za) m f°° dq f(dz)2-\ 



7T n 2 Jo {-q 2 ) 2 J U 2 IV UV° 

For the diagrams i?04 an d H12, the product of two electron propagators just outside 
the fourth order self-energy diagram behave as (1/q 2 ) 2 , which makes the convergence of the 
numerical integrals difficult in the small |g| region, even though the integrals are analytically 
free from the IR singularity in |g| after the mass and wave function renormalizations are 
carried out. In order to avoid this computational difficulty, we introduced an additional 
parameter y varying from zero to one to combine the original term and the renormalization 
term. All the numerator expressions are then proportional to at least q 2 and kills one of the 
electron propagators. 

The parametric representation of Group 3 is of the form 

3 a 2 (Za) m r°° dq f (dz)^ 1 

Ff~3. L / rro 775 ■ l 7 9) 



16 7T 7T 2 7o q 2 J U 2 V 3 

For instance, the diagram H 05 has the electron-line operator Fg'^: 

= 7 a (A + A + m) lp (p 3 + m)-f( A + m)<f(Ps + m)j a . (80) 

By using the hyperfine splitting projection operator, we get 



a 2 (Za) m r°° f {dz) x ^ 
Az/ G 3 = Ep - / dq 

7T 7T Z JO 



Fq _ F\ F 2 



V 3 UV 2 U 2 V 



(81) 



u 2 

The Group 4 diagrams Hi 7 , if 18 , and iJ 19 contain two non-overlapping second-order ra- 
diative corrections. Their sum is invariant under the covariant gauge transformation and 
free from the IR singularity due to the radiative photons. 

Let us consider the sum of two diagrams of Fig. |^. If the photon propagator is chosen 

as 
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FIG. 8. Sum of the second-order self-energy diagram and vertex diagram. 



(82) 



k 2 + it k 2 

the gauge-dependent part of the vertex diagram gives 

r A d A k i j -i 

J (27r) 4 j6— jk-\- fa — m fa— jk — m k 2 

= #f<$e* ^filA-rn «-V- (83) 
where we use the on- shell condition fp = m. The self-energy diagram with one external 
photon diagram is 

7 (27r) 4 fa — m fa+ fa — m k 2 



. , /A rf 4 fc 

^' (2^F A 



1 Y - Y 



(84) 



- Jk+ fa — m fa+ fa — m 

The second term, which is related to the mass renormalization constant proportional to the 
longitudinal photon polarization, vanishes when the integration over k is carried out with a 
proper regularization. Then the gauge dependent parts of fl55|) and ( |S4"D cancel each other 
and the sum is independent of particular choice of gauge. 

The numerical integration is performed for the integral combining three diagrams 
if 17, His and if 19 together so that cancellation of IR divergences occurs in the same re- 
gion of the Feynman parametric space. The result obtained for the zero mass radiative 
photon (A 2 = 0) is 

a 2 (Za) 

Au(H 17 ) + Au(H 18 ) + Au(H 19 ) = -0.478 03 (15)—^ -E F . (85) 

TT 
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This is in good agreement with the result calculated in the Fried- Yennie gauge [40|, in which 



P = -2 in (82) 



a 2 (Za) 

Au(H 17 ) + Au(H 18 ) + Av(H 19 ) = -0.477 89 (1) — -E F . 



7T 



(86) 



Note that Au(Hn), Az/(if 18 ), and Av(H 19 ) individually are gauge dependent, and their 
values are completely different between our results and those of Ref . (|(| . 



C. Problems Concerning Numerical Integration 

Let us now discuss some technical details of calculation of Au(H i) to Au(Hiq). After 
the ultraviolet divergences are renormalized, individual diagrams still suffer from severe 
infrared (IR) divergence, which is of the form A -1 / 2 , A being the photon rest mass measured 
in units of the electron mass. Of course, the sum over all diagrams of Fig. ^ is free from 
the IR divergence. This does not mean, however, that the sum can be integrated easily on 
a computer. This is because the IR finiteness results from cancellation of divergences for 
A —>■ from different parts of the integration domain. 

One way to deal with this problem is to evaluate individual integrals for several small 
values of A and extrapolate the sum of all terms to zero photon mass. Unfortunately, this 
approach creates integrals of order 10 3 for A 2 ~ 10~ 7 , while their sum is of order 1, making it 
very difficult to control the numerical accuracy of the result. Another way is to integrate, for 
A ^ 0, the sum of all terms, which enables us to avoid dealing directly with large numbers. 
This approach will also result in a better error estimate. The main practical difficulty is the 
large amount of computing time required. 

This problem can be somewhat alleviated if one evaluates each integral after subtracting 
its IR-divergent part, and then evaluates the sum S of the IR-subtraction terms of all 
diagrams. This method, which we have chosen, ensures that all integrals stay small (less 
than ~ 20) for any value of A. Thus far, we have evaluated them for several values of 
A 2 in the range of 10~ 3 to 10 -7 . The integration has been carried out numerically using 
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Table II: Vegas integration performed for individual diagrams of Fig. 2. S is the sum of IR divergent parts of diagrams 
Hqi through H w . X is the photon mass measured in units of the electron mass. 





A — 1U 


\2 _ 1 n-3. 5 
A — 1U 


A — 1U 


\2 _ 1 n-4.5 

A — 1U 


\2 _ 1 n-s 

A — 1U 


\2 _ m-5-5 

A — 1U 


\2 _ 1 n-6 

A — 1U 


\2 _ 1 n-7 

A — 1U 




-4.9551(4) 


-5.5548(4) 


-6.0481(4) 


-6.4490(4) 


-6.7714(5) 


-7.0304(4) 


-7.2360(4) 


-7.5298(4) 


H02 


0.0950(1) 


0.3601(2) 


0.6337(1) 


0.8976(2) 


1.1402(2) 


1.3565(2) 


1.5443(2) 


1.8398(2) 


H03 


-6.5527(4) 


-7.5536(5) 


-8.4213(5) 


-9.1585(6) 


-9.7758(5) 


-10.2877(4) 


-10.7089(4) 


-11.3314(4) 


H04 


1.6048(1) 


1.9811(2) 


2.3670(2) 


2.7574(3) 


3.1487(3) 


3.5389(3) 


3.9266(3) 


4.6937(3) 


H05 


4.0899(4) 


4.6116(5) 


5.0464(6) 


5.3997(5) 


5.6836(6) 


5.9083(5) 


6.0833(5) 


6.3199(6) 


J J 

-"06 


-3.4781(5) 


-3.8533(5) 


-4.2074(6) 


-4.5364(6) 


-4.8346(6) 


-5.1006(4) 


-5.3326(4) 


-5.7030(6) 


TT 


1.2588(3) 


1.2690(3) 


1.2426(3) 


1.1898(3) 


1.1176(4) 


1.0318(3) 


0.9361(3) 


0.7263(4) 


ZJ 

-"08 


-0.2167(0) 


-0.2980(1) 


-0.3677(1) 


-0.4255(1) 


-0.4721(1) 


-0.5092(1) 


-0.5380(1) 


-0.5775(1) 


Hog 


0.0059(2) 


-0.3047(2) 


-0.5929(2) 


-0.8486(2) 


-1.0682(2) 


-1.2528(2) 


-1.4061(2) 


-1.6333(2) 


Hio 


7.7184(3) 


9.1640(4) 


10.4746(5) 


11.6210(5) 


12.5985(6) 


13.4159(4) 


14.0900(4) 


15.0840(6) 


Hn 


-1.8437(2) 


-2.4663(3) 


-3.0333(3) 


-3.5309(3) 


-3.9559(3) 


-4.3127(3) 


-4.6078(3) 


-5.0457(4) 


Hu 


4.4780(3) 


5.4636(4) 


6.3444(5) 


7.1092(5) 


7.7567(5) 


8.2968(4) 


8.7423(4) 


9.4015(6) 


H13 


-3.2369(3) 


-3.6631(4) 


-4.0220(5) 


-4.3157(6) 


-4.5505(5) 


-4.7339(4) 


-4.8732(4) 


-5.0485(5) 


Hu 


-3.6981(2) 


-4.3597(2) 


-4.9672(3) 


-5.5084(3) 


-5.9802(3) 


-6.3839(3) 


-6.7246(3) 


-7.2437(4) 


H15 


-1.9611(2) 


-2.5301(2) 


-3.0774(2) 


-3.5843(2) 


-4.0403(3) 


-4.4414(2) 


-4.7882(3) 


-5.3325(4) 


Hi6 


-1.5242(1) 


-1.7852(1) 


-2.0058(2) 


-2.1863(2) 


-2.3307(2) 


-2.4442(2) 


-2.5325(1) 


-2.6519(4) 


Hn, Hig, H w 


-OA 


1780(1) 


S 


8.2098(17) 


9.4538(21) 


10.5296(27) 


11.4341(31) 


12.1829(36) 


12.7885(55) 


13.2507(104) 


13.8516(186) 


TOTAL 


-0.4840(20) 


-0.5436(25) 


-0.5828(31) 


-0.6128(35) 


-0.6295(40) 


-0.6381(57) 


-0.6526(105) 


-0.6585(186) 



FIG. 9. The graph for the coefficient y of E F a 2 (Za)/ir obtained by VEGAS versus x = A 2 . 
The solid line is the chi-square fit for y = a>o + a\x + ai x 2 , where ao,ai, and ai are determined 
from the data listed in Table I. 
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TABLE III. The chi-square fitting for the coefficient y of Epa 2 (Za) /it versus x = A 2 where 
y = ao + a\x 2 + • • • and A is the photon mass measured in units of the electron mass. 



photon mass A 2 




y = ao + clix 


y = ao + a\x + a2X 2 


y = ao + a\x + a^x 1 + a$x 3 




«o 


-0.678 ±0.015 


-0.700 ±0.041 





10- 7 < A 2 < 1(T 5 


di 


0.88 ±0.29 


2.03 ±2.05 







a-2 




-13.77 ±24.06 






a Q 




-0.6764 ±0.0079 


-0.672 ±0.017 


1(T 7 < A 2 < 10- 3 


a 1 




0.73 ±0.15 


0.58 ±0.57 








1.99 ±0.63 


3.45 ±5.62 




03 






-4.39 ± 16.78 



the adaptive iterative Monte-Carlo integration routine VEGAS [37]. The result of each 
integration is summarized in Table II. The degree of difficulty of numerical integration for 
S increases rapidly as A decreases. This prevents us from going to smaller values of A at 
present. Even A 2 = 10 -7 is a struggle. Although evaluation of integrals up to A 2 = 10 -8 
is highly desirable, we have not attempted it thus far since it will require an extraordinary 
amount of computing time. 

The data in Table II shows that the contribution of Fig. 1(f) falls within errors on a 
straight line for 10~ 7 < A 2 < 10~ 5 . Thus the extrapolation to A = may be tried with a 
linear polynomial a + cqx, where x = A 1//2 |39|. The upper box of column 3 of Table [HI 
shows the best linear fit to the data in this range of A. 

We also list in the upper box of column 4 a fit to the same set of data in terms of a 
quadratic polynomial ao + a\x ± a2X 2 . Clearly the result is much less certain than the linear 
fit. This is because the increased flexibility of the quadratic polynomial is now responding 
not only to the nonlinearity of data but also to the noise of numerical integration. For this 
reason the fitting with a linear polynomial will be more appropriate for this data. 

On the other hand, if one tries to fit the entire data of Table II which show clear deviation 
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from linearity, the linear fit is no longer appropriate and one must use (at least) a quadratic 
polynomial (See Fig. 9). The best fit by a quadratic polynomial to the whole set of data 
of Table II is shown in the bottom box of column 4. The bottom box of column 5 shows 
that use of a cubic polynomial is not recommended to this data because it responds more 
to the noise of numerical integration than to the real signal. Based on these considerations 
we believe that the value of ao determined by the quadratic fit gives the best estimate of 
AKFig.l(f)): 

Az/(Fig.l(f)) = -0.676 4 (79) ^ Za) E F . (87) 



V. DISCUSSION 

The remaining uncertainty in the value of A^(Fig.l(f)) is still considerable. Nevertheless 
it is a factor 5 improvement over the preliminary value fllCf). Since we published the pre- 
liminary result, Eides et. al. ]4T)|j4"l"|j have completed their calculation and reported a more 



accurate value 



Ai/(Fig.l(f)) = -0.671 1 (j) a2 ( Za) E F 

7T 



-0.370 1 (4) kHz. (88) 



Our new result (|87f) is in close agreement with the result of Ref. ||41|| . Note that the calcula- 



tion of Ref. |^T| is carried out in the Fried- Yennie gauge while our work is carried out in the 
Feynman gauge. Considerably higher precision of (|88|) over (^7j) reflects the advantage of 
calculation in the Fried- Yennie gauge which renders each diagram IR- divergence-free when 
it is transformed by application of integration by part. This approach requires nontrivial 
amount of diagram-specific manipulation because of very complicated integrands involving 
many variables. In the Feynman gauge calculation, on the other hand, each diagram is IR 
divergent due to the retarded Coulomb-like photons, which makes individual integral cutoff 
dependent. The advantage of this approach is that one can apply a systematic computer 
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algebraic method, which minimizes the chance of making mistakes — an important consid- 
eration in such a complicated calculation. The agreement of fl3?D and (|88D confirms gauge 
independence of the result to the extent of numerical precision. 

The total contribution of order a 2 (Za) including the results ([TT]) — (p~5|) and fl3T|) is given 
in flT8"|). If we instead use fl8"g| ) for Fig. [TJ(f), the total contribution becomes 

a 2 (Za) 

Az/(Fig.l) = 0.773 2 (7) 1 ' E F 

= 0.426 4 (4) kHz. (89) 

If we add the a(Za) 2 correction to the previous theoretical prediction (|i~0[), we obtain 



4 463 302.69 (1.34) (0.21) (0.16) kHz Eq.(18), 
Az/(new theory) = (90) 

4 463 302.70 (1.34) (0.21) (0.16) kHz Eq.flU). 

Now that the complete a 2 (Za) correction has been evaluated, the major remaining source 
of theoretical uncertainty in the muonium hyperfine structure is, as is seen from (|4]), the 
numerically evaluated nonlogarithmic part of the aiZa) 2 correction term. In addition, 
although not listed in (|90|), the leading logarithmic corrections of order a 4 and a 3 (m/M) 
turn out to contribute to the hyperfine structure as much as the a 2 (Za) term, as was shown 



in our preliminary report H and also by Karshenboim [42]. The parts of these higher order 



terms evaluated thus far add up to —0.68 (6) kHz if some errors in || and |TIJ are corrected. 
(Details will be discussed in subsequent papers.) We cannot simply add these corrections to 
(p0|), however, because the previous evaluation of aiZa) 2 term |23| contains parts which are 
of higher order in Zot. We have evaluated the aiZa) 2 and a(Za) 3 terms separately in the 
NRQED formalism in order to avoid possible double counting and reduce the uncertainties 
due to these terms. The results will be reported in the subsequent papers. 
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APPENDIX A: PROJECTION OPERATORS FOR QED RENORMALIZATION 

CONSTANTS 

We present the projection operators for QED renormalization constants determined by 
the on-shell scheme. Our projection operators for QED renormalization constants are quite 
general and correct for any loop orders. We set the electron mass m — 1 through Appendices 
and |B|. 

The projection operator of vertex renormalization terms can be written as 

L= 1 -Tr[T ( 1 ° + l)] J (Al) 

where is the proper vertex diagram. Projection operators of mass renormalization term 
and wave function renormalization terms can be defined similarly. A proper self energy 
diagram of the 2n-th order has the parametric representation 

E(2n) = " (^)V - 2 ) !F / • ( A2 ) 

(See Ref. [|35| for the definition of U, V, F, etc.) The projection to the mass renormalization 
is given by 

6m=^Tr[E(p)^+l)]\^ =1 

_ fay f (dz) G F m 

Kir) J U 2 ^ U m V n - 1 -™' { ' 

where is defined by ( |A2|) . The projection operator of the wave function renormalization 
constant is slightly more complicated since it involves the derivative of the self energy with 
respect to the external momentum, namely p' 1 dTi(p)/dp fl . It can be easily realized, however, 
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in the Feynman parametric representation. Taking the derivative of the numerator with 
respect to p leads to the quantity 

E= A i F i- ( A4 ) 

electron line only 

where Fj is obtained by replacing ( / D i + 1) by jp in the electron-line operator F. The 
corresponding part of the wave function renormalization constant is determined as 

( -a\ n l r 1 

burner.) = -(-) (n - 2)!^^) ftfi+l)]^! J (d e ) __ r 

a\ n f (dz) G n ^ E m 



ttJ J u 2 ^u m v n - l - m 

The other part corresponding to the derivative of the denominator is similar to the mass 
renormalization 8m and is given by 

B,( W= (g)7^g ("-^l Gf " , (A6) 



7T 



where 



2 F dp" 



p 2 =1 



£ *4 • (A7) 

electron line only 



APPENDIX B: NUMERICAL CALCULATION OF THE a(Za) CORRECTION 

Let us begin by showing the calculation of the spanning photon diagram of Fig. |5](a). 
The Feynman parameters assigned to the electron-lines are Zi,z 2 , and z 3 and assigned to 
the radiative photon is z±. Following the notation in Ref. [[34|,[JIJ, we use the abbreviation 
for the sum of Feynman parameters 



2l2-.fi = zi + z 2 H V z n . (Bl) 

For this diagram, Z1234 = 1. The momentum flowing through each electron- line after inte- 
gration of photon loop momentum should be expressed by the external momenta, namely 
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the electron momentum p = (m, 0) and the momentum of the external potential q — (0, q). 
For the electron-line i, this may be written as 

Q'i = Aip + A iq q, (B2) 

where the "scalar" currents Ai and A iq are found to be 

At = A 2 = A 3 = 1 - z 123 U~\ A lq = A 3q = -z 2 U~\ 

A 2q = 1 + A lq , U = z l23A = 1 . (B3) 

Then the electron-line including all of numerical factors is obtained in the operator form 



3 

8 q 



^[( A + 1) Y (A + 1) i u (A + 1)] / , (B4) 



where the i-th electron line operator Di is 

Df = i / dm 2 ^- . (B5) 

Multiplying with the hyperfme splitting projection operator and taking the trace, we obtain 
the integrand as the FORM output PE||: 



/(<?> 4 A 2? gV- 2 + (-4 A\ A 2q + 8\-4 A 2(/ ) 

+5i 2 (8i l5 + 4 A 2q ){UV)-\ (B6) 

where 

B 12 = B 23 = B 31 = l, V = z 123 - z 123 At + z 4 X 2 + q 2 z 2 A 2q . (B7) 
The contribution to the hyperfme splitting energy is obtained as an integral of the form 

AEt = a(Za)E F ± j™ dq J ( ^^f(q), (B8) 

where 

(dz)i-4 = dz\dz 2 dz 3 dz/Jb{\ — z\ 234 ) , (B9) 
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and f(q) is given by flB6|). 

Next consider the vertex correction diagram of Fig. |5](b). Feynman parameters Z\ and 
z 2 are assigned to the electron lines and Z4 to the photon line, with the constraint Z124 = 1. 
The scalar currents are 

A 1 = A 2 = l- z l2 U'\ A 2q = -z 2 U~\ A lq = l + A 2q , U = z 124 = 1 (BIO) 
Then the electron-line is 

-f ^[(A + l)/(A+l)7^+ai)]/ A f z,dmlJ { ^^j- 2 . (Bll) 
The integrand is found to be 

f(q)= -4 A l3 i 2g gV- 2 + (-4 - 4 A 1 A 2q + 4 A x A lq + 8 A^V' 2 

+B 12 (~A)(UV)-\ (B12) 

where 

£12 = 1, V = z 12 - z 12 Ax + z A A 2 + cf 2 z 2 A 2g . (B13) 



As was described in Sec. |HI|, the charge renormalization and subtraction of the second 
order anomalous magnetic moment is carried out by subtracting f(q = 0). The hyperfine 
splitting contribution from two diagrams of the vertex type is 



AE 2 = a{Za)E F ^ J™ £ z,dm\ J { ^±[f{q) - /(0)]. (B14) 



If the m\ integral is performed, the denominator V~ 2 becomes V -1 . The V 1 term is UV 
divergent and, if it is combined with the charge renormalization term, it becomes — In V. 
In the self-energy diagram calculation, the regularization is implicitly performed and the 
resulting denominators are expressed by V~ n , n — 0, 1, 2, - • where V~° implies — In V. 

The last diagram is the self-energy insertion diagram of Fig. |5](c). Feynman parameter 
z 2 is assigned to the electron-line and z 4 to the photon line with the constraint z 24 = 1. The 
scalar currents are 
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A 2 = l-z 2 U-\ A 2q = A 2 , 17 = 324 = 1. (B15) 
Then the electron-line is 

I ^[(^+^ + l)7 M (A + l)7 l/ (i»+^ + l)]/^ i ^, (B16) 

where 

V = 22 - z 2 A 2 + 2 4 A 2 + q 2 z 2 A 2q . (B17) 
The integrand is found to be 

f(q) = (-16 + 8 A 2 - 4 A 2 q 2 )(- \n(V)) . (B18) 
The mass and wave function renormalization terms can be written as 

f R = G(1Q - 8 A 2 )q Vq- 1 + (-16 + 8 A 2 - 4 A 2 <f 2 )(- ln(\/ )) , (B19) 

where 

G = z 2 A 2 , V = z 2 - z 2 A 2 + z 4 X 2 . (B20) 

The hyperfine splitting contribution is given by 

AE, = a (Z a) E F ^ fj^-J y^V® - hi ■ (B21) 

Although this integral is analytically free from the IR singularity caused in the limit of the 
vanishing external photon momentum q, numerical integration is very difficult in the small 
q region, resulting in the poor convergence of the integral. To avoid this numerical difficulty 
we introduce another parameter y varying from zero to one and combine the term which is 
not proportional to q 2 in the original diagram and the corresponding renormalization term 
together. This leads to an integrand of the form 

f(q) = G(16 - 8 A 2 )q W- 1 - 4 A 2 q 2 (- \n(V)), (B22) 

where 
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Vy = z 2 - z 2 A 2 + z 4 \ 2 + q 2 y z 2 A 2q . (B23) 
The corresponding renormalization term is 

f R = G(16 - 8 A 2 )q 2 V ' 1 - 4 A 2 q\- hi(V )) . (B24) 
The resulting hyperfine splitting contribution is 

AE 3 = a(Z«)*4 f ^ /; d y J ^[f ffl - /«] . (B25) 
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